The zero-field superconducting piiase transition obscured by finite-size effects in thick 

YBa2Cu307 5 films 
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We report on the normal-superconducting phase transition in thick YBa2Cu307-5 films in zero 
magnetic field. We find significant finite-size effects at low currents even in our thickest films 
(d = 3200 A). Using data at higher currents, we can unambiguously find Tc and z, and show 
z = 2.1 ± 0.15, as expected for the three-dimensional XY model with diffusive dynamics. The 
crossover to two-dimensional behavior, seen by other researchers in thinner films {d < 500 A), 
obscures the three-dimensional transition in both zero field and the vortex-glass transition in field, 
leading to incorrect values of Tc (or Tg), v, and z. The finite-size effects, usually ignored in thick 
films, are an explanation for the wide range of critical exponents found in the literature. 
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PACS numbers: 74.40.-|-k, 74.25.Dw, 74.72.Bk 



Since the discovery of high-temperature superconduc- 
tors and the realization that their higher critical tempera- 
tures and smaller coherence lengths create an experimen- 
tally accessible critical regior^i, researchers have looked 
at these superconductors in an effort to determine the 
model that governs the phase transition. Fisher, Fisher 
and Huse^ codified the scaling approach to the normal- 
superconducting (N-S) phase transition and predicted 
the existence of an N-S phase transition in field, called 
the vortex-glass transition. This phase transition in field 
has been extensively studied using current- voltage (I—V) 
isotherms, and although a consensus has emerged that a 
vortex-glass transition exists, there is little consensus in 
the values of the critical exponents v and a^i^. More- 
over, some have claimed that scaling data collapse does 
not prove the existence of a phase transition^ and that 
screening can create a non-zero resistance^, destroying 
the transition. Recent work has questioned the existence 
of a phase transition, showing that data collapse alone 
is too flexible, and proposing a criterion to determine 
whether or not a phase transition has occurred^. 

In zero magnetic field the existence of an N-S phase 
transition is not in doubt. Very close to {\T — Tc\ < 
2 K)\ the transition is expected to obey the three- 
dimensional (3D) XY model, with ly « 0.67 and z — 2.0 
for diffusive dynamics^.. Although specific heat and pen- 
etration depth measurements have found mean-field val- 
ues of ~ 0.5^, others have fit specific heat and penetra- 
tion depth data using critical models with smaller resid- 
uals than for mean-field modelsiS, and recent thermal 
expansivity data is more consistent with 3D-XY scaling, 
(i/ « 0.67)ii. Transport measurements can determine 
both and z, but data in zero field is inconsistent: Re- 
searchers have found vortex-glass like exponents — 1.1, 
z = 8.3) in small fields (< 10 mT)iS, others finding 3D- 
XY-like exponents when extrapolating to zero field from 
higher fields"'^'^ and in crystals^'*. 

Researchers have shown that, in thin films {d < 500 
A), the fiuctuation dynamics can cross over from D — 3 
to D = 2, and that this crossover occurs at a well-defined 



10' io^J<^"^o^ 




FIG. 1: I - V curves for a 2100 A YBa2Cu307_5 film, 
with bridge dimensions 20 x 100 fiia^ , in zero magnetic field. 
Isotherms are separated by 60 mK. The dashed line indicates 
a slope of 1, or ohmic behavior. The error bars are smaller 
than the points. The inset is R{T) at 10 ^A. 



current density, J^jrt^i?^. In this work we present a sys- 
tematic study of Jmin in films of different thicknesses. 
We find that even in our thickest film {d = 3200 A) the 
crossover to I? = 2 obscures the phase transition, causing 
incorrect choices for v, and z. However, at currents 
greater than Jmin, we see behavior as predicted by scal- 
ing which gives reliable values for Tc and z. 

We examined the zero-field I — V curves of 
YBa2Cu307_5 (YBCO) films deposited via pulsed laser 
deposition onto SrTiOs (100) substrates. X-ray diffrac- 
tion verified that our films are of predominately c-axis 
orientation, and ac susceptibility measurements showed 
transition widths < 0.25 K. R{T) measurements (inset to 
Fig. ^ show Tc « 91.5 K and transition widths of about 
0.7 K. AFM and SEM images show featureless surfaces 
with a roughness of « 12 nm. Our films also have a high 
critical current (Jc(77 K) « 2 x lO^^A/m^). These films 
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FIG. 2: dLogE/dLogJ vs. / for the / - V curves from Fig.[T] 
Isotherms are separated by 60 mK. The conventional choice 
for Tc, 91.26 K, is clearly not a horizontal line. The opposite 
concavity criterion can be seen at higher currents (/ > 40 
nA) about 91.44 K. The inset shows the 91.26 K isotherm for 
three bridge widths on the same film: 20 /xm (solid line), 50 
(dashed line) and a 100 ^m (dotted fine), which do not 
agree as a function of /. 

are of similar or better quality than most YBCO films 
reported in the literature. 

Our films were photolithographically patterned into 4- 
probe bridges of varying widths (8 - 200 fim) and lengths 
(40 - 1000 /im) and etched with a dilute solution of phos- 
phoric acid without noticeable degradation of R{T). We 
surround our cryostat with /i-metal shields to reduce the 
ambient field to 2 x 10~^ T inside. To reduce external 
noise, the cryostat is placed inside a screen room with 
low-pass T filters at the screen room wall and low-pass tt 
filters at the cold end of the probe. 

Figure n shows I —V curves taken on a 2100 A thick 
film on a bridge of dimensions 20 x 100 /im^. Scaling 
predicts^ 

ve+'-''ii^x±{ie-'iT), (1) 

where D is the dimension, z is the dynamic critical ex- 
ponent, ^ is the coherence length, and x± are the scaling 
functions for above and below the transition tempera- 
ture Tc- Fluctuations are expected to have a typical size 
^ which diverges near Tc as ^ ~ \T/Tc — ll^"^, defining 
the static critical exponent v. 

Above Tc at low currents, the I~V curves are expected 
to be ohmic (represented in Fig. ^ as a dashed line with 
slope 1), whereas at higher currents the isotherms are 
expected to show non-linear, power law behavior (slope 
greater than 1). Exactly at Tc, the coherence length di- 
verges while the voltage remains finite, which is true only 
if y oc (for D = 3), i.e. a straight line on a 

log- log plot^'^. Conventionally, Tc is chosen as the first 
isotherm without an ohmic tail, the isotherm at 91.26 K 



in Fig. n Data at higher currents and voltages are typ- 
ically excluded from fits because it is assumed that the 
system is being driven out of thermal equilibrium. The 
thick solid line at 91.26 K is a fit to a power-law at lower 
voltages, and gives a dynamic exponent z — 5.5, simi- 
lar to exponents found elsewhere^^, but clearly not the 
expected z = 2. 

We have suggested^ that a better way to determine 
the critical isotherm is to examine the derivatives of 
log£' vs. logJ isotherms. On such a graph, the critical 
isotherm would be obvious as a horizontal line with in- 
tercept (z -I- l)/2, separating isotherms with positive and 
negative slope (corresponding to concave up and down 
in Fig. nj . Our opposite concavity criterion^ states that 
isotherms at equal temperatures away from Tc should 
show opposite concavity at the same current level. 

The derivative plot for the I — V curves in Fig. ^ is 
shown in Fig. |3 There is no isotherm which is horizontal 
over the entire range of currents, contrary to theoretical 
expectations. The opposite concavity criterion is also not 
satisfied, and isotherms below 91.44 K have unexpected 
behavior: they are concave down at higher currents be- 
fore displaying ohmic behavior at lower currents. If we 
consider only the higher currents (/ > 40/iA) , we can see 
behavior as predicted by scaling which also satisfies the 
opposite concavity criterion: the isotherm at 91.44 K is 
horizontal, lower isotherms are concave up, and higher 
isotherms are concave down. This allows an unambigu- 
ous choice for Tc, 91.44 K. If we fit the high-current data 
to a horizontal line, then z = 2.1 ± 0.15, which agrees 
with diffusive dynamics. Below we will justify analyzing 
only / > 40/iA, ignoring the low-current linear behavior 
in these I — V curves (the ohmic "tails" ) . 

To determine whether the ohmic tails are a bulk in- 
trinsic effect, we patterned bridges of different widths 
on the same film. The inset to Fig. |21 shows the 91.26 
K isotherm for three bridges on the same 2100 A film 
from Figs. Hand 121 Each bridge was measured simulta- 
neously to insure identical temperatures. It is clear that 
the isotherms do not agree as a function of /. In Figure 
O we plot dlogE/dlogJ as a function of J rather than 
/. All three bridges have similar behavior in J, show- 
ing that we are measuring a bulk effect as opposed to an 
edge effeclii^. It is also clear that each isotherm turns 
over towards ohmic behavior at a certain applied current 
density rather than current. This is significant, because 
at an applied current density J one probes fiuctuations of 
typical sizeSiSa Lj = {ckBT/^oJf^, where $o = h/2e is 
the magnetic fiux quantum and c is a constant expected 
to be of the same order as the YBCO anisotropy param- 
eter, 7 « 0.2. Thus, as J decreases, Lj increases and will 
eventually reach the thickness of the film. At this point 
a crossover to 2D behavior is expected, as the size of the 
fiuctuations is limited along the c-axis. Thus, for a film 
of thickness d, there is a minimum current density, such 
that smaller current densities probe 2D fluctuations: 

Jmm = cksT/^od^. (2) 
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FIG. 3: dhogE / dLogJ vs. J for three bridges of different 
widths on the same 2f00 A film (20 x fOO^im^, 50 x 250fim^ , 
100 X 500/im^). The crossover to ohmic behavior clearly de- 
pends on J. 

Because Jmin doBS not depend on the exponents v and z, 
this minimum current density apphes for both the vortex- 
glass transition and the transition in zero field. 

We examined the ohmic tails generated in seven films 
of different thicknesses to determine how Jmin varies 
as a function of thickness. We measured seven films 
with similar properties (Tc and ATc) which varied in 
thickness from 950 A to 3200 A^^. To choose I -V 
curves to compare between films, we have taken the 
isotherm which, from high-current data, most seems like 
Tc, i.e. horizontal on the dlogi^/dlogJ vs. J plot. For 
Jmin we have chosen a similar criterion as Ref.^^, when 
dlog-E/dlogJ = ff the ohmic tails are caused by 

finite-size effects, then we expect l/^/JmAn vs. d to be 
a line with slope ($o/cfcsT)-'^/^, which will give a value 
for the undetermined constant c. The temperature T of 
the different isotherms only varies from 91.4 K - 92.5 K, 
a total change of about 1%. 

The results are plotted in Figure 01 Each value for 
1 / V Jmin incorporates error in Imin, bridge width and 
thickness, leading to error bars of about ±22%. Nonethe- 
less, the trend is clear: as d increases, 

Jmin decreases. 

The solid line in Fig. ^is a weighted least-squares linear 
fit to the data with a reduced chi-squared of = 0.41^^. 
From the slope we determine c = 0.60 ± 0.17, the same 
order of magnitude as 7 ~ 0.2, as expected. 

These finite-size effects have been seen by other re- 
searchers in thinner films. Dekker et ali^ found z — 
2.2 ± 0.4 from high-current data in a 500 A thick YBCO 
film in zero field, and saw ohmic tails at low currents. 
Using an equation nearly identical to Eq.|21 Dekker et al. 
noted that the fiuctuation size along the c-axis saturated 
at « 470 A, as expected. 

Finite size effects have also been seen in single crystals, 
both in a field and in zero field. Yeh et aL-^^ found ohmic 
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FIG. 4: 1/V J,nin vs. thickness for eight different films. The 
solid line is a weighted least-squares linear fit to the data, with 
a reduced chi-squared = 0.41. The slope of the line gives 
c = 0.60 ± 0.17, of the same order as 7 ~ 0.2, as expected. 



deviations from their data collapse at low currents and 
attributed them to finite size effects. The length derived 
from Eq. |21 agreed well with the distance between twin 
boundaries, and they suggested these boundaries limited 
the size of the fluctuations. Although they found good 
agreement between theory and experiment, the devia- 
tions were determined after the data collapse. We have 
shown that I — V curves can be made to scale with dif- 
ferent choices of Tg, v, and a^, thus apparent agreement 
with scaling via a data collapse is not conclusive evi- 
dence that a phase transition occurs, or that one's choice 
of critical parameters is the correct one. 

Woltgens et alM found deviations from 3D-scaling 
which appeared as ohmic tails in films with d < 500 A, 
as compared to a 3000 A film. Woltgens et al. assume 
that the finite-size effects in thin films do not extend 
to the 3000 A film because the I — V curves for the 
3000 A film scale with typical vortex-glass exponents, 
despite the fact that Yeh et al. found finite size effects 
in crystals, where the distance between twin boundaries 
were « 2 /xm, nearly an order of magnitude thicker than 
Woltgens et al.^s thickest films. Moreover, a simple data 
collapse is not conclusive evidence that the 3000 A thick 
films are unaffected by the finite-size effects they see in 
thinner films. For T = 83 K, J^in ~ 1 x 10'5A/m^ and 
data below Jmin were included in their analysis of the 3D 
transition, and included in the scaling collapse, possibly 
affecting their choice of Tg. This indicates that assum- 
ing a data collapse a priori and analyzing deviations from 
this collapse is not the correct method to determine finite 
size effects. 

Because the crossover to D = 2 can affect the choice 
for Tc, it is an explanation for the wide range of critical 
exponents found in the literature. I — V curves are ex- 
pected to be ohmic at low currents for T > Tc (or Tg), 
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thus it is possible to confuse ohmic tails generated by 
finite-size effects with ohmic tails generated by the 3D 
phase transition. This changes the conventional choice 
for Tc (or Tg) (the first isotherm without an ohmic tail), 
and because the ohmic tails are used to determine v (as 
R cx {T/Tc — f at low currents), then values for v 

and z will also be affected. For example, Sawa et al?^ 
scaled I — V curves for films as thin as 180 A and as 
thick as 10,000 A in a 2 T applied field by systematically 
changing Tg, v and z, citing this as evidence for the need 
for an anisotropic 3D-XY model. We suggest that the 
crossover to I? = 2, occurring at different current densi- 
ties, required them to vary v and z in order to scale their 
data. 

Our results indicate that low-current ohmic tails are 
due to finite-size effects. Failure to account for this leads 
to significant underestimates of Tc (or Tg) and incorrect 
values of v and z. We show in the derivative plot that 
the high current data agree with the opposite concav- 
ity criterion and lead to unambiguous choices of Tf. and 
z. Because the source of the low-current ohmic tails is 
in question, this leaves only the data collapse to find 
but using data collapse to find the critical exponents is 
perilous^. We can collapse the data using T^ and z found 
from the high-current data, which yields v w 1.2, simi- 
lar to values found elsewhereiSiSS, but this value for v is 
clearly not 3D-XY, and other values of Tc, z/, and z can 
collapse the data also. It is also unclear how the ohmic 
tails affect data at higher currents, especially when T is 
far from Tc and the critical region is small. 

It is also interesting to note that as d increases, Imin 
does not decrease without limit. Because / = J{wd), 



Imin = {w/d){ckBT/^o), usiug Eq.[21 The smallest Imin 
can be is when w = d, or Imin ~ 1 x 10^^ A (T = 90 
K). Thus any applied current below 0.1 fiA will probe 
fluctuations limited by the thickness of the sample, inde- 
pendent of whether the sample is a thin film, thick film, 
or single crystal. 

In conclusion, we have looked at YBCO microbridges 
of various widths (8-200 /im) in seven films of different 
thicknesses (950-3200 A) whose zero-field I — V curves 
are consistent with low-current ohmic tails created by 
finite-size effects, even in the thickest films. In contrast, 
the behavior at currents greater than Jmin {I > 40/j,A in 
our film) agrees with the opposite concavity criterion as 
predicted by scaling, and gives the expected 3D-XY dy- 
namic exponent of z = 2.1 ± 0.15. Because finite-size ef- 
fects are usually ignored in thicker films, we suggest that 
the low-current ohmic tails thought to be the expected 
behavior for T > Tc are actually generated by finite-size 
effects at temperatures T > Tc and T < Tc- This effect 
will obscure the phase transition in all films, both in zero 
and non-zero magnetic field, leading to incorrect results 
for the critical exponents and temperatures. 
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